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SUMMARY 


Time histories are presented that illustrate the large motions 
which have been encountered in flight tests of some of the present-day- 
fighter airplanes. 

Results of some analog stvidies are discussed which indicate that 
variations in certain of the airplane stability derivatives could have 
an appreciable effect on these vmdesirable motions. 
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wing area, sq ft 
wing span, ft 
chord, ft 

mean aerodynamic chord, ft 

dynamic pressure, Ib/sq ft 

Mach number 

pressure altitude, ft 

sweep angle, deg 

aspect ratio, b^/S 

total aileron deflection, deg 

stabilizer deflection, deg 
mdder deflection, deg 
aerodynamic forces along X body axis, lb 
aerodynamic forces along Y body axis, lb 
aerodynamic forces along Z body axis, lb 
direction cosines 

aerodynamic rolling moments, Ib-ft 
aerodynamic pitching moments, Ib-ft 
aerodynamic yawing moments, Ib-ft 

linear velocities along X, Y, and Z body axes, ft/sec 
angular velocities about X, Y, and Z body axes, radians/sec 
disturbances in linear velocities 

angle of attack of X body axis to relative wind, deg 
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0 

0 





sideslip angle, deg 
pitch angle, deg 
bank angle, deg 
yaw angle, deg 

rolling-moment coefficient due to sideslip 

longitudinal static derivative, per deg 

damping- in-pitch derivative, per radian 
directional stability derivative, per deg 
damping- in-yaw derivative, per radian 
nondimensional pitch frequency squared 
nondlmensional yaw frequency squared 


Subscripts: 


0 

initial value 

max 

maximum value 

res0 

refers to resonance 


refers to resonance 


A dot over a symbol indicates 


in pitch 
in yaw 

differentiation with respect to time. 


RESULTS AND DISCUSSION 


Several of the fighter- type airplanes presently being flown have 
exhibited a strong coupling between their lateral and longitudinal modes 
of motion, particularly in rolling maneuvers. This coupling has led in 
some instances to very large, essentially uncontrollable airplane motions 
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in pitch and yaw. Two airplanes currently being flight tested by the 
National Advisory Committee for Aeronautics and which have experienced 
violent coupled motions are shown in figure 1. The airplane on the 
right is an unswept- wing airplane of aspect ratio 5*09 and that on the 
left is a swept-wing design of aspect ratio ^.^ 6 . Both of these 
airplanes have their mass concentrated primarily in the fuselage, as 

evidenced by the values of the Inertia parameter 


%- 


shown in the 


figure for the respective airplanes. This parameter approaches unity 
for designs which have extremely low inertia in roll relative to the 
inertia in pitch and yaw. For comparison, for fighters of World War H 
this parameter was of the order of 0.5 or 0.4. 


Results of a recent flight of the unswept-wing airplane are pre- 
sented in figure 2. The maneuver shown is an abrupt aileron roll from 
level flight at M = 1.05 at an altitude of 50,000 feet. The param- 
eters plotted are sideslip angle, angle of attack, rolling velocity, 
and the control deflections applied by the pilot. A sideslip angle of 
20° was obtained in this maneuver, and, in the l/2-second time interval 
between 4.5 seconds and 5.O seconds, the angle of attack changed from 
-15° to +19°. This change in a corresponds to a change in normal 
acceleration from roughly -6g to +7g. Results of a similar flight of 
the swept-wing airplane are presented in figure 5- The maneuver is an 
abrupt aileron roll from level flight at M = O.7O at an altitude of 
52,000 feet. The maximum sideslip angle encountered was in excess of 
-25°, and the angle of attack, although not recorded below -l6° because 
of the range of the measuring instrument, was estimated to have been 
larger thaii -25°. These large motions impose high loads on both the 
airplane and pilot, and hence are dangerous as well as undesirable. A 
point of interest is that the initial variations in a and p are 
different for the two airplanes discussed although the rolling velocity 
for both cases is negative. This difference Is attributable to the 
fact that the longitudinal principal axis of inertia is initially above 
the flight path for the swept-wing airplane whereas for the unswept-wing 
airplane this axis is initially below the flight path. 

The purpose of this paper is to bring to the attention of the 
autopilot designer the existence of these violent coupled motions, and 
to discuss the results of some analog studies which indicate that the 
severity of the pitch-yaw divergences which have been encountered might 
be alleviated by use of various types of automatic stabilization or 
regulation. These studies were carried out by NACA personnel at the 
Langley Aeronautical Laboratory and at the MCA High-Speed Flight 
Station. 
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In the theoretical approeujh to the problem it vas necessary to use 
nonlinear equations of motion to represent the airplane dynamic charac- 
teristics. These equations, referred to airplane body axes, are as 
follows: 


m [u + q (wo+ Aw) - ry] • Fx '*^*3 
m [v +r (uq+ A u) - p(wo+ Aw)] « Fy + Wmj 
m [w + gy - q (uo Au)] = F^ + Wn^ 

P ■ ^xz ~ ^y) ~ ^xz pq * 

ly q - (iz-^x) p»'-^lxz Ixe ■= M' 

izr-ixzP^- (iY-ix) pq-»-ixzq'' -Ixe^^eq-N' 

*3 ’fTijr-njq 
rhj* njp-ijr 

"3 ° *3P‘'^3P 


These general equations of airplane motion have been modified to include 
the gyroscopic moments due to the rotating engine. For purposes of com- 
parison the nonlinear terms, which normally are neglected in the analysis 
of small motions, are underlined. These terms will generally be most 
traportant for airplanes which have low inertia in roll relative to that 
in pitch and yaw, and which are capable of fairly high rolling velocities. 
Inclusion of these nonlinear coupling terms was necessary in this inves- 
tigation and the solution of these equations had to be obtained from 
some computing machine such as an analog computer. 

In order to determine if the large motions presented in figure 5 
for the swept-wing airplane could have been predicted from these non- 
linear equations of airplane motion, this particular flight run was 
simulated on an analog conputer. In this simulation it was assumed 
that the component of the airplane velocity along its X-axis was con- 
stant, which reduced the problem to five degrees of freedom. Estimates 
of the airplane's mass and aerodynamic characteristics were obtained 
from available wind-tunnel tests, flight data, and theory. Figure 
shows the comparisons between the flight record and the analog results. 
Flight results are shown as dashed lines and the calculated motions as 
solid curves. The agreement between the flight results and the analog 
calculations for the time interval shown indicates that the simulation 
was sufficiently accurate to predict the existence of the large motions 
encountered in flight. In order to obtain the agreement shown in this 
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figure ; it was necessary to make some adjustments in the preliminary 
estimates of the airplane stability derivatives. In particular, it 
was necessary to take into account the variation with angle of attack 
of the rolling-moment coefficient due to sideslip basis 

of this agreement, additional calculations were made for various amounts 
of aileron deflection for this particular flight condition so as to 
obtain an indication of the response characteristics of the airplane 
in aileron- induced rolls. The basic input was a ramp aileron deflection 
which was held deflected until the airplane rolled to a specific bank 
angle and the aileron was then returned to zero deflection. The rudder 
and elevator were assximed to be fixed in these maneuvers. The bank 
angles for which results were obtained were ±90° and +360°. 

Results typical of those gotten frbm these calculations for left 
rolls of 90° and 560° are presented in figure 5* The aileron deflec- 
tion in each of these cases was such as to give a maximum rolling 
velocity of -1.5 radians per second. The quantities tabulated in these 
and similar runs were the maximum rolling velocity, maximum sideslip, 
and maximum disturbance in angle of attack during the rolling maneuver. 
For cases in which the values of and ta. obtained when the aileron 

deflection was reduced to zero, which is the recovery phase of the 
maneuver, exceeded the values calculated in the maneuver, these values 
were also tabulated. A summary of the results obtained over the air- 
plane rolling-velocity range is given in figures 6 and 7 • Results for 
left rolls are presented in figure 6 and results for right rolls are 
presented in figure 7. The maximum changes in sideslip angle and angle 
of attack are plotted against the maximum rolling velocity obtained in 
the rolling maneuvers. The difference between the left and right rolls 
is due to the asymmetric engine gyroscopic moments which were mentioned 
earlier. The magnitudes of Ax and A3 in the 90° rolls are seen to 
be fairly small. The worse condition shown appears to be the 560° left 
roll and the subsequent discussion will be limited to this case. As 
the aileron deflection, and hence the maximum rolling velocity, is 
increased, there is a consistent increase in the maximum variations in 
sideslip angle and angle of attack for rolling velocities up to -2.6 
radians per second. The values of approximately 24° for A3 and l4° 
for ta were obtained for this roll rate. For p\irposes of comparison, 
the values of rolling velocity which produce resonance in pitch and yaw, 
based on a steady- rolling analysis (see refs. 1 and 2), are shown by 
the dashed vertical lines. Some preliminary analysis of the steady- 
rolling case has indicated that variations in certain of the airplane 
stability derivatives could have an appreciable effect on the magnitudes 
of a and 3 encountered in these rolling maneuvers. Figure 8 presents 
divergence boundaries for the swept-wing airplane. These boimdaries, the 
construction of which is discussed in references 1 and 2, are plotted ^ 
a function of the airplane pitch and yaw natural frequencies nondimension- 
alized to rolling velocity. The solid curve is for the basic airplane. 

It can be shown that the position in this plane of the pitch and yaw 



NACA RM L 55 E 251 ) 


7 


frequencies of a given airplane, nondimensionalized to rolling velocity, 
falls on a straight line similar to the one shown in the figure. The 
slope of this line is determined from the ratio of the square of the 
pitch and yaw natviral frequencies. For p = 1 radian per second the 

values of and (x>^ are equivalent to their dimensional counter- 

parts, and this point is spotted on the straight line in order to show 
the dimensional val\aes of the pitch and yaw frequencies of the swept- 
wing airplane. As the rolling velocity increases, the airplane frequency 

2 

parameters move down this line and for. a value of coy = O.7O cross the 

divergence hoimdary. For the swept-wlng airplane, this crossing of the 
bo\andary occurs for a rolling velocity of I.7 radians per second. The 

bovmdary is again crossed for w 1 which corresponds to a rolling 

velocity of 2.2 radians per second. For values of rolling velocity 
between I.7 and 2.2 radians per second, this simplified analysis pre- 
dicts a divergence, and it is these values which were indicated by the 
vertical lines in figure 6. 

It is apparent from figure 8 that this divergent condition can be 
eliminated in at least two ways. Variations in the airplane dimensional 
pitch or yaw frequencies can rotate this line to prevent its passing 
through the unstable region, or the boundaries can be moved with respect 
to the line which represents the swept-wing airplane to effect a similar 
result. These boundaries can be shifted by increasing the damping in 
pitch or yaw as indicated by the dashed boundaries for increased damping. 

For the condition shown, this line can be rotated into the stable 
region by increasing the yaw natioral frequency through increases in Cn^ 

or decreasing the pitch natural frequency through a decrease in 

However, for too large a variation in these parameters, the line will be 
rotated to a position for which it will pass through the lower branch 
of the divergence boundaries. It appears that a good rule of thimib 
would be to rotate the line to a position where the airplane pitch and 
yaw frequencies are approximately equal. 

Results will be presented in the subsequent figures which indicate 
the effect of variations in the static derivatives Cn^ and Cn^ and 

in the daniping derivatives Cn^. and Cases will be discussed for 

increases in Cno which rotate the line clockwise from its basic orien- 
P 

tation to the stable region and then to a position for which it passes 
through the pitch divergence boundary for a rolling velocity approxi- 
mately equal to the pitch natural frequency. For a variation 

will be discussed which rotates this frequency locus from its basic 
position to a position in the stable region of this plot. Cases will 
be discussed for Cnj. and Cm^ which shift the boundaries in the direc- 
tion indicated by the dashed lines (see fig. 8), and hence put the basic 
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frequency locus in a stable region. In the subsequent figures, values 
of rolling velocity are plotted which correspond to the resonant values 
predicted by the sin5)lified steady- rolling analysis. 


effect of varying the directional stability parameter Cnp is 

presented in figure 9. This derivative could be varied artificially by 
deflecting the rudder proportional to sideslip angle. The curve for 
Cn = 0.001 corresponds to the airplane without automatic controls. 

As Cnp is increased to 0.002, the maximum variations in sideslip and 

angle of attack computed for these rolling maneuvers are markedly reduced. 
It appears, however, that an Increase to 0.00l^ results in a condition 
during recovery which is about as bad as that which exists for 
Cnp = 0.001. For this value of Cnp, the frequency locus has been 

rotated to a position such that it passes through the lower bomdary 
for a rolling velocity approximately equal to the pitch frequency, and 
there is seen to be a large positive peak in Aa which exceeds the 
negative values calculated for the maneuver. It is apparent from fig- 
ure 9, however, that the large variations in angle of attack and side- 
slip experienced by this airplane in rolling manettvers are largely due 
to a basic deficiency in directional stability and, if the existing 
value of Cnp of 0.001 could be increased to a somewhat higher level 

by artificial means or by redesign of the vertical tail, the airplane 
characteristics would be improved. 


Results are presented in figure 10 for variations in the longitud- 
inal static derivative The curve for = -0.006 corresponds 

to the basic airplane. As Cidq^ is reduced to -0.005, for which case 

the frequency locus is in the stable region for all values of p, the 
a- and p-motions are seen to be substantially reduced. 


Results are shown in figure 11 to show the effect of the damplng- 
in-yaw derivative Cnj.. This derivative could be varied by use of a 

"yaw damper," in which the rudder is deflected proportional to yawing 
velocity. The value of Cnr estimated for the swept-wing airplane for 

the flight condition being discussed is -0.10. On the basis of the 
estimated effectiveness of the edrplane rudder, it appears that an 
increase in Cnj. to -1.0 is approximately of the order of magnitude 

which could be realized from a yaw damper for this flight condition. 

As Cnr is assumed to be varied from its basic value of -0.10 to -1.0, 

there is seen to be an improvement in the airplane characteristics in 
these aileron- induced rolls, although the inq>rovement is less than that 
calculated for siid 
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Results are presented in figure 12 for variations in the damping- 
in-pitch derivative Cmg. This derivative could be varied by a pitch 
dan 5 )er which deflects the stabilizer proportional to pitching velocity. 
As Cniq is ass^nned to be varied from its basic value of -5*5 to -55-0, 

the maximum value of ^ is reduced from 2lt-° to 6° and the angle-of- 
attack variation is reduced from ll<-° to 6°. This value of Cmq is 
approximately in the range which could be realized from a pitch damper 
for this particular airplane and flight condition. 


CONCLUDING REMARKS 


It should be en 5 )hasized that these preliminary calculations were 
made for only one flight condition and for one type of maneuver, and 
therefore no attempt should be made to generalize the results presented 
to other conditions of altitude and Mach number and to other maneuvers. 
However, the trends indicated by these results are believed to be 
qualitatively representative of those which would be obtained through 
use of the automatic controls discussed. 

The primary conclusions which can be drawn from this paper do, how- 
ever indicate that (l) certain of the fighter- type airplanes which are 
being flown today are highly susceptible to divergences in angle of 
attack and sideslip during rolling maneuvers and (2) people in the auto- 
matic control field may be in a position to assist in the solution of 
this problem by equipping airplanes with autopilots which will reduce 
the severity of these undesirable motions. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May 11, 1955- 
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CHARACTERISTICS OF AIRPLANES 
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Figure 2 
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RESPONSE OF SWEPT-WING AIRPLANE TO AILERON DEFLECTION 
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Figure 5 


COMPARISON OF FLIGHT AND ANALOG COMPUTER 
RESULTS FOR SWEPT-WING AIRPLANE 
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Figure 4 
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RESPONSES OF SWEPT AIRPLANE IN 360® AND 
90® LEFT ROLLS 
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Figure 5 


SUMMARY OF RESPONSES OF SWEPT-WING AIRPLANE 
TO AILERON DEFLECTION 
M*0.70; Hp* 32,000 FT 
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Figure 6 
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SUMMARY OF RESPONSES OF SWEPT-WING AIRPLANE 
TO AILERON DEFLECTION 

M = 0.70; Hp = 32,000 FT 

A^, 0E6 

90 

360 


Aa 

OR 

-A/9, 

DEG 


I I I ^ ^ U U I I 

0 .4 .8 1.2 1.6 2.0 2.4 2.8 3.2 

Pmax* radians/sec 

Figure 7 



DIVERGENCE BOUNDARIES - STEADY ROLLING 
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Figure 8 







